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Abstract 

Dark energy and the accelerated expansion of the universe have 
been the direct predictions of the distant supernovae la observations 
which are also supported, indirectly, by the observations of the CMB 
anisotropics, gravitational lensing and the studies of galaxy clusters. 
Today these results are accommodated in what has become the con- 
cordance cosmology: a universe with flat spatial sections t = constant 
with about 70% of its energy in the form of Einstein's cosmological 
constant A and about 25% in the form of dark matter (made of per- 
haps weakly interacting massive particles). Though the composition 
is weird, the theory has shown remarkable successes at many fronts. 

However, we find that as more and more supernovae la are ob- 
served, more accurately and towards higher redshift, the probability 
that the data are well explained by the cosmological models decreases 
alarmingly, finally ruling out the concordance model at more than 
95% confidence level. This raises doubts against the 'standard can- 
dle '-hypothesis of the supernovae la and their use to constrain the 
cosmological models. We need a better understanding of the entire 
SN la phenomenon in order to extract cosmological consequences from 
them. 
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1. Introduction 

The history of cosmology has probably never witnessed such upheavals as 
evidenced in the past few years, which are primarily due to the observations 
of distant supernovae (SNe) of type la. Many colourful exotic models have 
been speculated, which claim to provide satisfactory explanation to these 
observation. It is generally believed that the distant SNe la observations 
predict an accelerating expansion of the universe powered by some hypothet- 
ical source with negative pressure generally termed as 'dark energy'. The 
simplest and the most favoured candidate of dark energy is a positive cos- 
mological constant A, which is though plagued with the horrible fine tuning 
problems. This has led a number of cosmologists to resort to scalar field 
models of evolving dark energy, which can produce negative pressure for a 
potential energy-dominated field and cause the scale factor to accelerate at 
late times by violating the strong energy condition. While the scalar field 
models enjoy considerable popularity, they have not helped us to understand 
the nature of dark energy at any deeper level. 

Though the idea of the accelerating expansion and dark energy in the 
framework of Einstein's theory is a prediction of the first generation of SNe 
la data, however as more and more accurate data get accumulated, the fit 
to the cosmological models worsens successively and the recent observations, 
taken at their face value, seem to rule out all the cosmologies at fairly high 
confidence levels! This will be shown in the present paper. It may be noted 
that in our goodness-of-fit analysis, we include only those observations which, 
unlike the SNLS observations, have already included the intrinsic dispersion 
of the SN absolute magnitude in their error bars. The SNLS data are not 
suitable for a goodness-of-fit analysis, as we shall see later. 

It should be noted that the precise measurements of the temperature 
anisotropies of the CMB made by the WMAP experiment [T], which appear 
to offer the most promising determination of the cosmological parameters, 
are often quoted as providing a direct evidence for an accelerating universe, 
which is though not quite correct. The standard interpretation of the WMAP 
constraints may be misleading, as it relies on the assumption of the power 
law spectrum. Blanchard et al. [2] have shown that the CDM Einstein- 
de Sitter (EdS) universe is quite consistent with the WMAP data if the 
primordial spectrum is not scale-free. Hence, taken on their face values, the 
only apparent prediction made by the WMAP observations is a flat geometry, 
and the decelerating models like the EdS also explain these observations 
successfully [3J. 

2. m — z relation in Robertson- Walker cosmologies 

As the measured quantities of the SNe are the magnitude (m) and the redshift 
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(z), let us derive, in brief, the usual m — z relation for SNe la in the framework 
of Einstein's theory, which we shall use for our analysis. The derivation 
assumes the simplest model of the universe based on the Robertson- Walker 
(R-W) metric, representing a homogeneous and isotropic spacetime. Let us 
assume that the observer at r = and t = to receives light emitted at t = t\ 
from a SN of absolute luminosity L located at a coordinate distance r x . The 
measured (apparent) luminosity I of the SN is defined by I = L/Airdl, where 
g?l, the luminosity distance, is given by 

d L = {l + z)S n, (1) 

where z = 5 (to) / S(ti) — 1 is the cosmological redshift of the SN, with S being 
the R-W scale factor. (Incidentally the luminosities I and L are expressed in 
terms of the K-corrected mag nitudes m and M as I = 10" 2m/5 x 2.52 x 10" 5 
erg cm -2 s" 1 and L = 10~ 2M / 5 x 3.02 x 10 35 erg s" 1 .) When written in 
terms of the magnitudes m and M, the above luminosity-redshift relation 
gets transformed into the magnitude-redshift relation: 

m{z; M, Sli) = M + 5 log{H d L (z; (2) 

where M. = M — 5 log 10 H + constant. The value of the constant depends 
on the chosen units in which d\, and Ho are measured. For example, if g?l is 
measured in Mpc and Ho in km s _1 Mpc -1 , then this constant comes out as 
~25. The value of T\ appearing in flTJ can be calculated by integrating the 
R-W metric for the SN: 

sin (A So 7^) ) » when k = 1 

iSoW^l' when k = (3) 

. sinh (i So W&)), when k = -1. 

The Hubble parameter H(z), appearing in these equations is provided by the 
Einstein field equations. If the different sources which populate the universe 
do not interact with each other and each of them is represented by an equation 
of state Ui = Pi/ pi (which can be a function of time in general), the Friedmann 
equation then yields 



H\z) = Hi 



E^exp UJ q Z l \^ ) dz> \ - n k (1 + z) 



(4) 



where Qi are, as usual, the present day energy densities of the different source 
components in units of the critical density 3Hl/8nG and fife = k/S^H^ 
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(i denoting non-relativistic matter (m), radiation, (r), cosmological constant 
(A), quintessence ((f)) etc.). The present value of the scale factor So, appearing 
in equations (Ofl [3]) which measures the present curvature of spacetime, can 
now be calculated from 

We note from equations (J31 [3J) that the coordinate distance r 1 , and hence d^, 
are sensitive to fij for the distant SNe only. For the nearby SNe [z << 1), 
equation (j2J) reduces to 



m(z) = M. + 5 log z, 



(6) 



which can be used to measure M. by using low-redshift supernovae-measurements 
(that are far enough into the Hubble flow so that their peculiar velocities do 
not contribute significantly to their redshifts). 

Now for given M., fli and u>i, these equations can provide the predicted 
value of m(z) at any given z. We compare this value with the corresponding 
observed magnitude m and compute x 2 horn 
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N 

E 



m(zj] M., Qi, ujj) — m Qj , 
cr m .„ , 



(7) 



where the quantity «7 m . is the uncertainty in the observed magnitude m j 
of the j-th SN. It may be noted that sometimes the zero-point absolute mag- 
nitudes are set arbitrarily in different data sets. While fitting the combined 
data set this situation is handled successfully by the constant M. appearing 
in equation (J7J), which now plays the role of the normalization constant and 
simply gets modified suitably. In this case however it does not represent the 
usual 'Hubble constant-free absolute magnitude' but differs from the latter 
by an unknown constant (which is though not needed for the cosmological 
results). The constant A4 also takes care of the cases where the data are 
given in terms of the distance modulus \i = m(z) — M, instead of m. Equa- 
tion d?D can also be used in this case for fitting the data by using \x a in place 
of m . 

Sometimes in the data we are also provided with independent uncertain- 
ties on some another variable (say, y). In this case the equation for \ 2 g e ^s 
modified as 
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{m(zj] M., Vti 
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where a y j is the uncertainty in the observed variable y corresponding to the 



The key point about the SNe la data-fitting is that the absolute lumi- 
nosities M of all the SNe, distant or nearby, are regarded same (standard 
cand /e-hypothesis) . Hence so is the constant Ai, as it has only one extra 
parameter H which certainly does not differ from SN to SN. Thus there are 
two ways of the actual data fitting: (i) estimate M. by using equation (jSJ) 
from the low-redshift SNe, and use this value in equation (j2J) to estimate fi; 
from the high-redshift data alone; (ii) use low-, as well as, high-redshift data 
simultaneously to evaluate all the parameters from equation (TSJ) by keeping 
Ai as a free parameter. Obviously the second method gives a better fitting 
(which we have used throughout this paper). 

It is obvious from equation ((7|) that if the model represents the data cor- 
rectly, then the difference between the predicted magnitude and the observed 
one at each data point should be roughly the same size as the measurement 
uncertainties and each data point will contribute to x 2 roughly one, giv- 
ing the sum roughly equal to the number of data points N (more correctly 
N— number of fitted parameters = number of degrees of freedom 'dof'). If 
X 2 is large, then the fit is bad. However we must quantify our judgment and 
decision about the goodness- of- fit, in the absence of which, the estimated pa- 
rameters of the model (and their estimated uncertainties) have no meaning 
at all. An independent assessment of the goodness-of-fit of the data to the 
model is given in terms of the x 2 -P r °bability: if the fitted model provides a 
typical value of x 2 as i at n dof, this probability is given by 



P(x, n) gives the probability that a model which does fit the data at n dof, 
would give a value of x 2 as large as x or larger. If P is very small, the model 
is ruled out. For example, if we get a x 2 — 20 at 5 dof for some model, then 
the hypothesis that the model describes the data genuinely is unlikely, as the 
probability P(20, 5) = 0.0012 is very small. 

3. Fitting cosmological models to different available SNe la data 

sets 

Now we have developed enough infrastructure to analyze the observations. 
We start our analysis with the data from Perlmutter et al. [I], which is 
one of the important data sets of the first generation of SN la cosmology 
programs [U [3] . We particularly focus on the sample of 54 SNe from their 
'primary fit' C. We have shown the fit-results in Table 1. The ACDM, as 
well as the models with the constant have a good fit. For example the 



j-th SN. 




(9) 
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concordance model (flat ACDM) has a x 2 — 57.7 at 52 dof with the goodness- 
of-fit probability P = 27.3%, representing a good fit. The Einstein de Sitter 
(EdS) model (O m = 1, A = 0), which used to be the favourite model before 
the SNe la observations, has a bad fit with P = 0.06%. In order to take note 
of the historj0, we have also shown the fit to the Bondi-Gold-Hoyle steady 
state model (for which m = M. + 51og[z{l + z}]) - the first model which 
predicted an accelerating universe. 

The existing data points coming from a wide range of different observa- 
tions were compiled by Tonry et al. [6]. With many new important additions 
towards higher redshifts, a refined sample of 194 SNe was presented by Bar- 
ris et al. [7]. However, the data we are going to consider next, is the 'gold 
sample' of Riess et al. [S], which is a more reliable set of data with re- 
duced calibration errors arising from the systematics. It contains 143 points 
from the previously published data, plus 14 new points with z > 1 discov- 
ered with the Hubble Space Telescope (HST). While compiling this sample, 
various points from the previously published data were discarded where the 
classification of the supernova was not certain or the photometry was incom- 
plete. The fit-results of this data show that although the fits to the ACDM 
and the constant u^-models are reasonable, they are deteriorated consider- 
ably; the probabilities P have reduced to less than half of the corresponding 
probabilities obtained in the case of the Perlmutter et al' data. Earlier, Roy 
Choudhury and Padmanabhan [U] have also shown that the Riess et al' 'gold 
sample' is inconsistent with a flat cosmology at 90% confidence level. 



2 A number of theories, e.g., brane cosmology, phantom cosmology, Cardassian cos- 
mology, Chaplygin gas cosmology, have been proposed as alternatives to the standard 
paradigm. However either they reduce to the standard ACDM cosmology in the present 
phase of evolution or have even worse fit than the standard ACDM cosmology . So, they 
are not included in the fit. 
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Models 




or Vt^ u><p 


M 


x 2 


dof 


P 




54 SNe from Perlmutter et al. (1999) 










ACDM (flat) 


0.28 ±0.08 


1 - fim "I 


23.94 ± 0.05 


57.7 


52 


0.273 


ACDM (n.c.) 


0.79 ± 0.47 


1.40 ±0.65 -1 


23.91 ±0.06 


56.9 


51 


0.266 


cons (flat) 


0.48 ±0.15 


l-n m -2.10 ±1.83 


23.91 ± 0.08 


57.2 


51 


0.257 


EdS 


1 





24.21 ±0.03 


92.9 


53 


0.0006 


Steady State 


2 





23.78 ± 0.03 


75.8 


53 


0.022 




Gold sample 


of 157 SNe from Riess et al. (2004) 








ACDM (flat) 


0.31 ±0.04 


i-n m -l 


43.34 ± 0.03 


177.1 


155 


0.108 


ACDM (n.c.) 


0.46 ±0.10 


0.98 ±0.19 -1 


43.32 ± 0.03 


175.0 


154 


0.118 


cons cj0 (flat) 


0.49 ± 0.06 


l-n m -2.33 ±1.07 


43.30 ± 0.04 


173.7 


154 


0.132 


EdS 


1 





43.58 ± 0.02 


324.7 


156 


io- 13 


Steady State 


2 





43.15 ±0.02 


318.3 


156 


10" 13 


164 SNe from Gold + ESSENCE (Krisciunas et al. 2005) 








ACDM (flat) 


0.30 ± 0.04 


i-n m -i 


43.34 ± 0.03 


190.3 


162 


0.064 


ACDM (n.c.) 


0.48 ±0.10 


1.04 ±0.18 -1 


43.32 ± 0.03 


187.2 


161 


0.077 


cons cj^ (flat) 


0.50 ±0.06 


l-n m -2.69 ±1.30 


43.29 ± 0.04 


185.5 


161 


0.091 


EdS 


1 





43.60 ± 0.02 


348.2 


163 


io- 15 


Steady State 


2 





43.15 ±0.02 


331.1 


163 


io- 13 


168 SNe from Gold + ESSENCE + 4 SNe from Clocchiatti et al. (2005) 






ACDM (flat) 


0.29 ±0.04 


1 - "I 


43.35 ± 0.03 


200.8 


166 


0.034 


ACDM (n.c.) 


0.51 ±0.09 


1.12 ±0.16 -1 


43.32 ± 0.03 


195.8 


165 


0.051 


cons (jjfj, (flat) 


0.50 ± 0.04 


l-n m -3.18 ±1.46 


43.28 ± 0.04 


193.4 


165 


0.065 


EdS 


1 





43.61 ±0.02 


367.5 


167 


io- 17 


Steady State 


2 





43.16 ±0.02 


338.3 


167 


io- 13 


Gold sample + 4 SNe from Clocchiatti et al. (2005) 








ACDM (flat) 


0.30 ± 0.04 


1 ^~^m 


43.35 ± 0.03 


188.0 




U.UOo 


ACDM (n.c.) 


0.49 ±0.10 


1.07 ±0.17 -1 


43.32 ± 0.03 


184.2 


158 


0.075 


cons (jJtj, (flat) 


0.50 ±0.05 


l-n m -2.95 ±1.41 


43.29 ± 0.04 


182.1 


158 


0.092 


EdS 


1 





43.60 ± 0.02 


345.0 


160 


10~ 15 


Steady State 


2 





43.16 ±0.02 


325.6 


160 


io- 13 



Table 1. Fits of different cosmologies to available data sets: some models have 
been constrained by the requirement of a flat space (f2 m + = 1), whereas the 
rest have no constraint (n.c). 

We now consider the first results of the ESSENCE project (TO] (made 
public in August 2005) under which 9 SNe with redshift in the range 0.5 - 
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0.8 were discovered jointly with HST and Cerro Tololo 4-m telescope. In 
order to minimize the systematic errors, all the ground-based photometry 
was obtained with the same telescope and instrument. We consider 7 SNe of 
this project which have unambiguous redshift and definite classification, and 
add them to the 'gold sample' resulting in a reliable sample of 164 SNe. The 
Table 1 shows that the fits to different cosmologies have further worsened 
considerably and do not represent a good fit, in any case. Increasing the 
number of fitted parameters (for example, in the models with a constant uj^) 
improves the fit marginally only. 

Next we consider the recent discovery of 5 SNe at redshift z ~ 0.5 by 
the High-z Supernova Search Team [11] (results made public in October 
2005). We consider 4 SNe from this sample for which distances estimated 
from the MLCS2k2 (Multi-colour Light Curve Shape) method are available, 
so that we can include them in the previous sample of 'gold + ESSENCE' 
which also use the MLCS2k2 method to determine the distance moduli. The 
fit-results of the resulting sample of 168 SNe are very disappointing. The 
quality of the fits to different models has deteriorated to such an extent that 
the concordance model can be rejected at 96.6 % confidence level! This is an 
alarming situation. Other models have marginally similar fit and increasing 
the number of fitted parameters does not help significantly. Models with 
variable cj^(t) do not help either. For example, if we consider lo^{z) = ujq + 
uji z/(l + z) with u , uj\ as constants, we obtain Q m = 0.42, = 0.42, 
ujq = —4.95 and ui = 2.83 as the best-fitting solution with x 2 — 193.07 at 
163 dof and P = 5.4% (in fact, the model is very degenerate in this case and 
the parameters wander around near the minimum \ 2 m a flat valley of some 
complicated topology). 

It has been mentioned that the ESSENCE data are affected by the selec- 
tion bias which can affect their cosmological use [TU]. However, we do not find 
any such effect from the fitting, as is clear from the Table. We notice that 
the parameters and their uncertainties estimated by using the ESSENCE 
sample are absolutely consistent with those from the use of gold sample or 
the gold+Clocchiatti et al.'s sample; and we do not find any smoking gun 
pointing out that this sample is completely different from the others. How- 
ever, even if we exclude the ESSENCE sample from the fit, the concordance 
model is still ruled out at more than 94% confidence level. 

Finally we consider the recently published (made public in October 2005) 
first year-data of the planned five-year SuperNova Legacy Survey (SNLS) 
|12j . However, the SNLS data have been analyzed in a different way than 
the other data sets we considered earlier and it does not make sense to add 
them for a joint analysis. Hence we limit ourselves to commenting on the 
way SNLS data have been analysed. 
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In the SNLS, the authors claim to have achieved high precision from 
improved statistics and better control of systematics by using the multi-band 
rolling search technique and a single imaging instrument to observe the same 
fields. Their data set includes 71 high redshift SNe la in the redshift range 
0.2 - 1 from the SNLS, together with 44 low redshift SNe la compiled from 
the literature but analyzed in the same manner as the high-z sample. As the 
SNLS data have been analyzed differently, the fitting procedure followed by 
the authors is also different. In order to calculate \ 2 i they use 

<j + a int I ' 

where a- mi , is the (unknown) intrinsic dispersion of the SN absolute magnitude 
which, unlike the other data setj§, is not included in the a^ ] rather it has 
been used as an adjustable free parameter to obtain x 2 /dof = 1. We shall 
return to this issue later for our comments. First we want to verify if we 
can reproduce the results of Astier et al. from equation (jiup by using their 
calculated \x a (given in columns 7 of their Tables 8 and 9) instead of using 
their 'stretch' and 'color' parameters, which do not seem necessary once we 
have Ho. We find that by fixing cr int = 0.13, we get fl m = 1 — Qa = 0.26, 
M = 43.16 with x 2 /dof = 1.00; and fi m = 0.31, fi A = 0.81, M = 43.15 with 
X 2 /dof = 1.01. This is exactly what Astier et al have obtained. 

The intrinsic dispersion in the absolute magnitude of SN la (combined 
with dust extinction of the host galaxy) can be estimated only statistically 
(unlike the photometric error, which can be estimated from the photometry 
of the individual SN la). Unfortunately we do not have a reliable way to 
estimate dust extinction or pure intrinsic dispersion of SNe la separately. 
However, the introduction of <j- int in equation (11 Op is justified only when 
we use independent measurement uncertainties (Tintj on the parameter (as 
we have mentioned earlier in equation (jSJ)), instead of using it as a free 
parameter. The latter case is just equivalent to increasing the error bars 
suitably in order to have a desired fit. In this way one can fit any model 
to the data. For example, the EdS model can also have a similar fit by 
considering cr int = 0.258: M = 43.46 with xVdof = 1.00. This shows that 
the approach does not have any predictive power. One could choose the 
variable uj^{z) such that it gives a lower x 2 /dof with the same a\ nt = 0.13. 
Also, for a similar o"i nt , one can obtain a reasonable value of the reduced x 2 
for another model. For example, increasing cr irit only to cr int = 0.16, the model 

3 It should be noted that all the other observations, we have considered before, already 
include the intrinsic dispersion of the SN absolute magnitude in their error bars a mo or 
cr jL , o , which has been estimated by reasonable methods 



x 2 = E 

7 = 1 
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Q m = 0, = gives x 2 /dof =0.86. However, one cannot physically test 
the viability of the model so that one can take any decision. This happens 
because the present approach (which simply assumes, rather than tests, that 
the data have a good fit to the model) prohibits an independent assessment 
of the goodness-of-fit-probability P, in the absence of which the estimated 
parameters do not have any significance. All one can do, with the present 
approach, is that one can compare the goodness-of-fit of different models. For 
example, with cri nt = 0.13, the EdS model has a worse fit than the ACDM 
model, giving x 2 /dof = 2.7. In this context, it would be remarkable that 
the results of the SNLS data are consistent with a flat geometry whereas the 
best-fit values of the other observations give a larger value of the fitotai- 

4. Conclusion 

Supernovae la observations have profoundly changed cosmology by predict- 
ing an accelerated rate of cosmic expansion, and thus a repulsive dark energy 
component - an issue which is regarded as an almost mature science now. 
However, as more and more accurate data get accumulated, thanks to the 
remarkable progress made in various types of astrophysical and cosmological 
observations in recent years, they do not seem to fit any cosmology. The 
recent observations, taken at their face values, seem to rule out all the cos- 
mologies at fairly high confidence levels. Though these probabilities may 
not be regarded sufficient to rule out the models completely, however, they 
are high enough to point out towards the alarming trend of the recent data: 
as you add newer data to the older samples, the goodness-of-fit-probability 
from the resulting samples successively decreases. Though the fit improves 
in some cases if we do not stick to the concordance model, however, this is 
inconsistent with the anisotropy measurements of CMB which predict a flat 
space. 

The recently made SNLS observations by the Supernova Legacy Survey 
[T2] are analyzed in a different way. While the other observations estimate 
the intrinsic scatter in the absolute magnitude cr^t from the nearby data, 
the SNLS estimates it from all the observed data. For this purpose, a- m t is 
considered as a free parameter to be estimated from all the data by requiring 
that it gives a good fit. Obviously, this data is not suitable for a goodness- 
of-fit analysis. It must be noted that our result (that the recent observations 
seem to rule out all the cosmologies at fairly high confidence levels) is deduced 
from those observations only which, unlike the SNLS data, have already 
included a- mt in their error bars. 

Assuming that the standard big bang cosmology is correct, the present 
situation is pointing out towards some flaws in our understanding of the SN 
la phenomenon and towards the futility of the use of SNe la in order to 
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constrain cosmological models. We need better understanding of the entire 
SN la phenomenon in order to test the empirical calibrations that are so 
confidently extrapolated at high redshifts. Similar conclusions have also been 
drawn by Clocchiatti et al. [H] from a smaller sample of data. However, this 
is more evident from the present analysis of a bigger sample of data. This 
view is also supported by Middleditch [T3] who argues that SNe la seem to be 
affected by some systematic effects which alone, without invoking any dark 
energy, could make them too faint for their redshifts. It is argued that it 
may be impossible to get a clean sample of SNe la which are free from this 
kind of effects [13] . 

This reminds us of a similar story of the m — z test for galaxies in 1970s - it 
was agreed upon finally that uncertain evolutionary effects tend to dominate 
at high redshifts. Though most studies confirm that the luminosity prop- 
erties of SNe la at different redshift and environments are similar [U [15], 
however, there are other theoretical studies which have found variations in- 
dicating evolutionary effects [16] . It has been shown by Drell et al. p2] that 
the peak luminosities estimated for individual SNe la by two different meth- 
ods are not entirely consistent with one another at high redshifts, z ~ 0.5. 
If evolution was entirely absent, the differences between them should not de- 
pend on redshift. They further showed that the three luminosity estimators 
in practice (the multicolor light curve shape method, the template fitting 
method, and the stretch factor method) reduce the dispersion of distance 
moduli about best fit models at low redshift, but they do not at high red- 
shift, indicating that the SNe have evolved with redshift. This view is also 
corroborated by some recent observations which may go against the 'standard 
candle'-hypothesis of SNe of type la, we mention the following two: 

1. SNLS-03D3bb, which is a recently observed high redshift (z = 0.244) 
type la SN with extreme unusual features and no obvious analogue at 
low redshifts. It does not obey the usual lightcurve shape-luminosity 
relationship for SNe la that allows them to be calibrated as standard 
candles [18] . 

2. Observations of two supernova remnants DEM L238 and DEM L249 
made with the Chandra and XMM-Newton X-ray satellites in the large 
Magellanic cloud [H] may also be mentioned. While the presence of 
Fe-rich gas at the centres of these objects suggests that they are rem- 
nants of type la supernova explosions, the standard model of type la 
supernova remnants cannot explain the presence of relatively dense su- 
pernova ejecta with long ionization timescales. 

We may recall that in the 1970s-80s astronomers tried to fit number 
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counts of radio sources to cosmological models and improved the fit by adding 
multi-parameter evolutionary functions like luminosity evolution, number 
density evolution, etc. This gave interesting possibilities of evolution though 
had no predictive value for the model. One can do similar exercise to im- 
prove the fit of the SNe la, particularly by keeping the dark energy equation 
of state parameter uj^(z) completely free as a function of z. 

Note added in the proof: It may be mentioned that a more recently 
made observation [20] claims improvements in the data and provides better 
fits to the cosmological models. However, this sample does not include the 
ESSENCE data [10] in the analysis, and hence our analysis is still meaningful. 
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